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ABSTRACT 

This paper describes an investigation of the dynamic behavior of 
a pin ended elastic column, subjected to half-sine pulse loading applied 
with small eccentricity. The column is replaced by a lumped parameter 
mathematical model, and the equations for the model are solved with a 
high speed digital computer. The failure critericn used is a limiting 
value of extreme fiber strain. The minimum loads which cause failure 
are found as a function of load duration for columns having the slender- 
ness ratios 50, 100, and 150. It is shown that an elastic column can 
Support rapidly applied dynamic loads greatly in excess of the Euler 
load. As the duration of the load pulse is decreased, the lateral 
deflection at failure becomes progressively smaller and the effects 


of axial inertia become increasingly significant. 
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1. Introduction 

The problems of the dynamic behavior of structures are ones which, 
like so many others, have only recently been attacked with any sort of 
vigor by the engineering profession. Whether this is due to the previous 
unimportance of the problems or to the lack of tools for their solution 
is primarily of historical interest; the fact remains that the problems 
are of importance today, and tools are now available for the solutions. 

One of these problems of structural dynamics which has aroused 
considerable interest in recent years is that of a dynamically loaded 
column, This seemingly simple structural member, whose behavior under 
conditions of static loading was predicted by Euler in 1757, becomes a 
quite complex system when the loading is applied dynamically. 

Attacks on this problem have been concentrated thus far on solu- 
tions for two types of loading - constant velocity loading of one end 
of the column, and impact loading. Hoff [1,2] : has treated the case 
of an elastic column, initially curved in the shape of a half sine wave, 
subjected to constant velocity loading such as that encountered during 
compression tests in commercial testing machines. He has shown that 
rapidly loaded slender columns with small initial deflections will 
support loads greatly in excess of the Euler load. Chawla les has 
extended this work to include the case of inelastic columns. Sevin| 4 | 
has confirmed Hoff's results, while retaining the effects of axial 
inertia (which were not considered by Hoff), and one of his conclusions 
is that 


»..SO long as the column remains elastic, axial inertia effects 
are of negligible importance in so far as the gross behavior of 


aNoMbers in brackets refer to bibliography. 
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conventional structural columns is concerned regardless of the 
initial deflected shape, end fixity, or type of axial loading. 


Gerard and Becker [5] have studied the impact loading case using the 
unloading strain wave produced upon failure of a tension specimen, and 
have concluded that a column may momentarily withstand any magnitude of 
compressive stress, and that the buckling may occur over a small portion 
of the length of the column, rather than the entire length. 

The literature also contains several treatments [6,7,8] of the 
problem of the stability of a pin ended column subjected to an axial force 
of the form Py = Py + A sinwt. It has been shown that, for certain 
values of w , the maximum compressive force se + A may become much higher 
than the Euler load without causing instability. On the other hand, it 
is also found that instability may exist when PY is a tensile force, 
provided that A and » have the proper values. 

Konig and Taub [9] have treated the case of a pin ended column with 
an initial half sine deflection, subject to a suddenly applied force of 
constant magnitude and variable duration. Their investigation shows 
that perfectly elastic columns can support loads in excess of the Euler 
load when the duration of the load is short. 

Other than in this last reference, the problem of the prediction of 
the load carrying capacity of a dynamically loaded column subjected to 
an externally applied force pulse of specified shape and duration seems 
to have been neglected. This is the problem which is considered in the 
present investigation, 

The column is assumed to be perfectly elastic, and is initially 
straight, rather than having some initial curved shape. It is of con- 


stant cross section, has constant physical properties, and is free of 





any damping. The effects of rotary inertia and shear strains are neglected, 
but axial inertia effects are retained, as well as non-linear axial strain 
components due to bending. 

The loading imposed on the column is a half sine pulse of force, 
applied with an arbitrarily chosen eccentricity, in an axial direction. 
The column is considered to have hinged-hinged end conditions, with the 
loaded end free to translate in the longitudinal direction. The unloaded 
end has an eccentric fixed-pin connection which allows only rotation. 

Since the concept of stability or buckling of the column seems to 
lose its meaning when applied to columns subjected to dynamic loading, 
some other criterion of failure must be used. In this study, an arbitrari- 
ly selected value of the extreme fiber strain is used to define failure. 

The problem is formulated and solved, not in terms of the real column, 
Ege in terms of a lumped parameter type of mathematical model. A set of 
equations is developed for this model and is then solved in a high speed 


digital computer. 


2. The Mathematical Model and D2velopment of the Equations for the System. 


In order to study the dynamic behavior of the column, the real column 
is replaced by a lumped parameter model consisting of a series of hinged 
rods, with point masses at the hinged joints. A general section of the 
mathematical model is shown in Fig. 1. The complete model consists of 
n increments, each of length £ = L/n, and n + 1 point masses, each having 
amass of m s Ah (except those at either end, whose masses are 2 £/2). 
The model is initially perfectly straight, with an eccentric fixed pin 
connection at mass n+ 1. The external load is applied eccentrically to 


mass number one, which is free to rotate and to move in the x direction, 





but is restrained from motion in the y direction. The eccentricity 
of the fixed pin connection at mass number n + 1 is the same as that 


with which the loading is applied. 





Figure 1. A section of the mathematical model of the column. 


The masses are connected by perfectly elastic, massless rods, 
which are hinged to each other at the point masses, as shown. The rods 
have an axial stiffness of K = EA/£% , where E is Young's modulus for the 
column material, and A is the area of the cross section of the column, 
The rods are considered flexurally rigid. 

At each of the point masses, with the exceptions of those at either 
end, is a perfectly elastic, massless spiral spring, with a spring con- 
stant C = EI/@ for either direction of rotation. (I is the rectangular 
moment of inertia of the cross section of the column) 

Referring now to the free body diagrams shown in Fig. 2 and Fig. 3, 


the following equations may be written: 








Fig. 2 Free body diagram of one mass of the mathematical 
model of the column. 
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Fig. 3 Free body diagram of one of the length increments 
(rods) of the mathematical model of the column. 





Fig. 4 The geometry of one of the Jength increments of the 
mathematical model of the column. 





arte) dt —an(P) seme (1) 
du /dty = (Vir Vin)/ mi (2) 
Ros + Vi sinQ = KARL; (3) 
AM, = M:-Mi-, = PR avi ~V,2+4u)) (4) 
In these equations, AJ, = £ - Xi (see Fig. 4), AU = W- Uy 


a; wut Ue 5 M: ae (8; ois > and t is real time, 


Solving eqs. (3) and (4) for P- and V- yields 


a AM; sin®@i + KAL: Z+Auz) (5) 
: L ad De 

a KALZ:Avi - AMicos Si 6 
Me L- AL: Z 


In order to put these into dimensionless form, both sides of eqs. 
(5) and (6) are now divided by F = Ee /t1* , where P. is the Euler load 


for the column. The numerator and denominator of both eqs. are divided 


by 2 , while noting that’ 
K/F = (EA/e) (ler) = (EAM) -("P'/ea) = MO/L 


where = £is the slenderness ratio cf the length increment, giving 


ee ee 
 * 1 — AL; /0 
and 
ce ‘iG 2 (ALi/9) (Avi /2) a (Al i eg) cosO: 


ee 
=e a - “ar 
——— @ —'* oP 


> , 0, 





The following approximations are now made; 
sinO; ~ ©; 
cosQ, vy 1-0:/2 
1 -AL/e ~ 4 
AV |p ~ Or 
These approximations are in error by one per cent or less for the 
values of Q@;, ALi, AUi, and AU; anticipated. 
Substituting in the equations for A; and By , and using the 


notation AM = AM: /FL , leads to 
Ac= AMidx + nia? (ALi/p)(14 AUvi/g) 


and 
B= rip (S4/p)e. - ami (1 - 9/2) 
Now, reference to the geometry of Fig. 4 shows that 
Lo = (L+ Auj)jcos%& + Aji sinO; 
or 
Ab =h-f' =L -(L+ AuJcos©; —- Avi sin ©, 
whence 


ALi ly - { — lest + Aui/y \cos G; — (Av:/p) sin & 


If the same approximations as were used previously are now substituted, 
this becomes 

ALi/y = Cize)(Auifg) - (Aui/y + Of/2) 
We further assume that the product (O%]2)(Au /2£) is negligible 
compared to the sum Ai /y =e Oi /2 , so that 

ALifg =~ - ( Aui fy Oi /2) 

This approximation is now substituted into the expressions for Aj and 
BL, giving 


Av= Am Or - mrp (AUi/y + Si7/2)(1 + AU¥i/s) 
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and 
Bi = -rtg? (Suc/g + e%/2) -—AM: (1 - %/2) 
These equations may be further reduced by assuming that 
(ay |AUife| << 4 
(b) |AMioc| << | npr ( Auily + Si /2) | 
(c) S/R << 4 
The physical significance of assumptions (a) and (c) is clear. 
The assumption made in (b) is equivalent to neglecting the contribu- 
tion of the y direction force component to the total force transmitted 
by the rod, i.e., omitting V, sin 6. in eq. (3). 
Based on the above, 


A, = -mpr (Atif + 9/2) 7) 


and 
Bs = Oc. A: -—- AYN: (8) 
Equations (7) and (8) are the dimensionless expressions for the 
components of the force transmitted by the ith length increment in the 
x and y directions, and are used in this form in the numerical solution. 
Using the notation AP; = P: —P;4; and AA; = AP, /e , equation 
(1) may be written as 
déui Jat = AA: /(mi/f) 
But, Mi ae = AL and F= EA | 0* , so that 
drui/dt®? = AALE/yf ro 
However, E / ee= Ss” » where S is the velocity of travel of an 


elastic wave in the column material, so that 


Pula = ah (* fem) = Bh(i(Ynp) © 


- 


- =. ¢ 4—— / i @ win 
 ialieeinen —_ <= - —— a -* | 
lc Terme io = fee 1 een 

— a 2 ee = af 


Se «e.'! 





and 


These equations may be further reduced by assuming that 
(a) 
(b) 
(c) 
The physical significance of assumptions (a) and (c) is clear. 
The assumption made in (b) is equivalent to neglecting the contribu- 
tion of the y direction force component to the total force transmitted 
by the rod. 
Based on the above, 
(7) 
and 
(8) 
Equations (7) and (8) are the dimensionless expressions components 
of the force transmitted by the ith length increment for the in the x 
and y directions, and are used in this form in the numerical solution. 
Using the notation and » equation 


(1) may be written as 


But, and , so that 


However, » where is the velocity of travel of an elastic 


wave in the column material, so that 


(9) 





For a beam with hinged ends the natural period of first mode vibra- 
tion ,7, , may be found from the expression 
2. Sou 4 Y2 
a =a b Jul?) 
Solving this for TU gives 
aed 20 np / 1s 
which leads to 
S*/p? z Ait 0° |? 7? 
Substitution of this expression in (9) yields 
Pui /gee a AML (4L |r [2) 
Or, since J is not a function of time, 
d(MYL) [yy = AAt (4% / ae) 
Dimensionless time is now defined as T = t/T, , so that 


(dt)* = T (dT). Therefore, 


LOM ype = AA (AM), 1eb4n — ca0) 


By the same reasoning as has just been applied to eq. (1), eq. (2) 

may be reduced to 
d*("Yk) [yr aN Gna , 14 en (11) 

Equations (10) and (11) are the equations for the dimensionless 
accelerations in the x and y directions respectively, for l1<i<£n, 
“which are used in the numerical solution. For the first mass, the 
acceleration in the x direction is twice the value given by eq. (10), 
due to this being a half mass. 

In order to "solve" the acceleration equations for the n + 1 masses, 
it is necessary to specify the forcing function, boundary conditions, 


and initial conditions. These are specified as follows: 





(a) for Of T#1/G6  , the external load applied (eccentrically) 
to mass number one is given by 


A a= F, sin(tt B T), where Pp = ayy ee 
ForT > 1/8 , A, = 0. Specifying the applied force in 
this manner allows a selection of both the amplitude and 
duration of the pulse, The bending moment at mass one is 
given by M, = P, e, while that at mass n+ lis PP, e. 
(b) forOf TS a; Vi/L=O 4 Uan/pzo , Un+i/p = O 
(¢) at T= 0, for l 2S aien + 2: 
Vi fy = 0 : Uilg = 0 ; dk (MY/2) 47 =O, d (“Phir = 0) 

The choice of the eccentricity with which the force is applied is 
made arbitrarily. The value used, given in dimensionless form, is e/r = 
0.05, where e is the actual eccentricity, and r is the radius of gyration 
of the column cross section, 

Another arbitrary choice which is made is that of a failure criter- 
ion for the column. For this investigation, the column is said to have 
failed when the extreme fiber strain, € 4» reaches a value of 0.01. 
This represents a stress of 300,000 psi in steel, which is admittedly 
high, but which is also an attainable yield point for certain alloy 
steels. It is felt that the use of this high value sets an upper limit 
for columns fabricated from presently available material. 

In order to compute the value of the extreme fiber strain, € ate it 
is necessary to find both the centroidal axis strain €,. and the strain 
caused by the application of bending moments, since C ax =lE.)+)e,| - 


The centroidal axis strain in any length increment of the model is 


calculated from 


Ect = ALi /y = — (Aui/g is Si /2) 
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The strain in the same length increment, caused by the bending moments, 
may be computed as 

ern = Ke Marg.) /Et 
where Marg. = (Mi + Mit) /2 
and c is the distance from the centroidal axis to the extreme fiber. 
Recalling that M; = FR Wi , we may say that 


Sain = (Fic/rer)(Mi + Mi-1) 


which reduces to 


CL (c/enea) (Mi + Mi-1) 
However, 2= rA , so that 


Et 


(¢/arrp)( Mi + Mi-i) 
It is also necessary to make a choice, at this point, of the value to 
be used for c/r. We know that for a thin walled, hollow cylindrical» 
cross section, c/r = (2)2 , while for a solid cylindrical cross section, 
c/r = 2. The value used in this investigation is c/r = 1.5. 
Substituting this value in the equation above gives 
Ey. = (0-75 /p*o) (Me + Yu) 


from which 


Ezmax = | AWi/y + B/2 | “+ | ozs (Yn + Mi-1)| (12) 





3. The Numerical Solution 

The equations of the system are solved by a numerical method of 
integration, utilizing Fortran programming and a Control Data Corpora- 
tion 1604 high speed digital computer. A simplified block diagram of 
the basic program is shown in Fig. 5, and a complete program, together 
with the program notation, is given in Appendix III. 


Referring now to Fig. 5, it will be seen that S, A, ,n, x, , 
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Fig. 5. Block Diagram of Basic Computer Program 







Main Program Subroutine Sigma 
sere see 
Read Compute 
‘ Vi 
5 Aio,1, 067 » 5 dO) LAT ¢! d{ M/s 
Set Initial 
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Compute 
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Compute Compute 
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Compute Compute 
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Compute Compute 
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Return 





and & are given as input data to the program. ( 8 is the ratio G/T, , 
where T, is the arbitrarily selected dimensionless time at which the 
program is to terminate). For the basic program shown in Fig. 5, the 
value of A, is that which corresponds to P, /P, = 1.0. 

At some time, T, , when the values of ce Ova A d’C/2)/aT*, 
ACYL) /AT , A (MEL) /AT 3 (iat) and (vi/2) have been found, 
the program jumps to the subroutine in order to compute these same 
quantities for time T, = T, + AT, where AT is the time increment used. 


As a first estimate, the accelerations are assumed to be constant during 


“seis dC 4/2) /AT 








the time increment, AT; that is, d*(*Y/2)/T? + 
1 


Based on this assumption, the velocity may be calculated as 
Ui = Ui + » (ui 4 
ACD AAT! = AMYM/dT| + ats d(HYQ/aT* |, 
(Oily the development for the x-direction quantities is shown here; the 


development is precisely parallel for the y-direction ) Having the velc- 


“| 


city, we now compute the displacement at T, as 





(ui/e)|. = (uso |. * (aT/2)fac“l/a7], + d(4/py AT 


which will be recognized as a trapezoidal integration method, 

With the first estimates of the displacements at T, now known, the 
forces A; and B; may be calculated from eqs. (7) and (8) (with the excep- 
tion, of course, of A, , which is specified). Eqs. (10) and (11) are 
then utilized to find the accelerations at T, . With these values in 
hand, the program now recomputes the velocities and displacements at T,, 


using the trapezoidal integration scheme: 


d (uilt)/aT |, = d (MJ) /dT | + (aT /2)) easy fat? |, + au) far, | 
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and 
(ur /O),, = GeO, + (aT/2)| cus /a7 ly + d(ui/e) fat, 

At this point, there is a return to the main program, and the Cima, 
are calculated from eq. (12). These Ci(may) values are then compared 
to the failure strain, 0.01. If any of the Ci(mayy 18 equal to or greater 
than 0.01, a printout of the quantities shown (in Fig. 5) is made and 
the program stops. If none is equal to or greater than the failure strain, 
the program compares the elapsed time with the specified terminate time, 
If the elapsed time is equal to or greater than the terminate time, the 
amplitude of the applied load is doubled, the initial conditions are 
reset, and the program commences again with the new value of applied 
load. If the elapsed time has not yet reached the terminate time, the 
computation is continued with the original value of A, . 

In this manner, a rough set of failure values of P, /P. is found. 
In order to refine these, a program such as the one given in Appendix III 
is used. It is basically the same as the one diagrammed in Fig. 5, with 
the following exceptions: 

(a) the input value of A, is that which was found to give 
failure using the basic program 

(b) the program does not stop after the printout of failure 
values. Instead, the A, value is decreased by a certain 
amount and a jump is made back to the point at which the 
initial conditions are set. The run continues until a non- 
failure value of A, is found. 

(c) the program stops after the time comparison, when T> 1/8, 


instead of continuing with a doubled value of A, . 
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Using this second program, any desired degree of refinement in 
the failure values of P, /P. may be attained. 

The program shows a marked sensitivity to the time increment, AT, 
used in the integration scheme. For 82 5, the time increment has been 
computed from 


Tt /2 
aay = n*P (0.48) 
which gives an equal number of time intervals during the force pulse, 


for any 8= 5. For &< 5, it was found that 


= Te 
AT = ae y. 


which corresponds in real time to At = SGeeris)., gives good accuracy (as 
judged from a work - energy comparison which is described in the follow- 
ing section). 

The computer program also shows some sensitivity to the number of 
length increments n used in the model. An investigation of this problem 
for values of n between 5 and 40 reveals that n = 20 gives good results 


without requiring excessive computer time, 





4, Discussion of Results 

Solutions for the failure value of P, /P, have been obtained for 
three values of S (column slenderness ratio), and for a range of @ from 
0.1 to 80. These results are summarized in Table 1 and are shown 
graphically in Fig. 6. 

It should be noted that the values of P, /P. given in the table and 
plotted in Fig. 6 are values midway between the minimum P, /P, 's giving 
failure and the maximum non-failure values. Adding the tolerance, © , 
to the tabulated P, /P, values gives the minimum failure P, /P.'s 


found, while subtracting 6 gives the maximum non-failure values. 
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It is possible to make some predictions regarding the values of 
P, /P_ which should obtain at both the high and low ends of the 8 range, 
and also, to predict the maximum values of P, [Pe . First consider the 
case for S = 0.1, which is approaching the static loading case. If the 


failure criterion of rem = 0.01 is applied to the secant formula for 


x 


eccentrically loaded columns, the following is obtained: 
- 2 P24: 
OC marx = 4 oe ==, Sec ($3 rig Sy ) | 


which becomes, on multiplication of both sides by A/P. 


mee = [t+ fabs) 


However, 
Omar A _ sy Fl = P sy 
=) 
so that 


a 
con ($) = [1 + $8 82¢(8fRS 9") ] 
Upon substitution of the selected values of e/r = 0.05 and c/r = 1.5, 
this becomes 
! 
E max (Sy = £/4 + Q.075 sec(E§-B2*) | 
This equation, when solved by a trial and error method, yields the follow- 
ing values of P, /Pe , compared here with the computer solution values 


for 6 = 0.1, for the three slenderness ratios: 





Slenderness Secant Computer 
ratio Formula Solution 
Value Value 
S P, /Pe P, /Pe 
50 0.943 0.925 
100 0.989 1.01 
150 0.996 1.03 


Another value of P, /P_ which may be predicted analytically is the 
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maximum load which the column can support without lateral deflection. 

If the column is loaded rapidly enough so that the amount of lateral 
deflection at failure is insignificant, then the column may be treated 
as an eccentrically loaded compression member, and the limiting load may 


be computed from 


C - _Tmax - {8/Pe)max 
ae St eeo AE (1+ 8S 
* Pz re 


which gives 
ea 
G = S a 
(Pi / Pe. ) avai = Emax (=) Cee a! 
The solution of this equation for S = 50, 100, and 150, — = OF 01e 


and ec/r* = 0.075 is compared below with the maximum values of Pires 


taken from the curves of Fig. 6. 


Slenderness Upper Limit Computer 
ratio S (P, /Pe ) max (P, /Pe ) max 
50 | 2.36 1.97 
100 9.42 78 
150 21.20 16.625 


In each case, the computed value falls below the limiting value. 
This is due to the fact that bending action is allowed in the computer 
solution, whereas the theoretical solution assumes that only axial com- 
pression takes place, and also, to the fact that the computer solution is 
based on dynamic, rather than static loading. 

The limiting value of P, /Pe which should be found at the high end 
of the 8 range may be predicted from a consideration of the travel and 
reflection of an elastic compression wave in the column. For a wave 
having the shape of a half sine pulse and a duration of T;/2, the 


following will be noted: 





(a) the wave front will travel from the loaded end of the 
the column to the pinned end in a time L/s. 


(b) at the pinned end, the original pulse will be reflected 
as a compression pulse traveling back toward the loaded 
end. This pulse will add to the original. 

(c) the compression pulse thus reflected will travel back to 
the loaded end of the column where it will again be re- 
flected, this time as a tension pulse which subtracts from 
the sum of the original and the first reflection. 

From this analysis it may be seen that if the loading is sufficient- 
ly rapid that the first tension pulse reflected from the loaded end does 
not arrive at the pinned end before the peak value of the original com- 
pressive pulse, then the maximum axial force at the pinned end will be 
twice the peak value of the applied load. 

Thus it is readily seen that, for Tr | + ae Pays (or B= G/ 3 > 
$/2%), the maximum value of P, /P. which the column can support will 
be one-half the (P, /P- )max value previously predicted, The applicable 


values of @ and -~(P, Pe ) axete given below for the three slender- 


ness ratios. 


50 7.95 1,18 
100 15.90 4.71 
150 23.85 10.60 


These limiting values of P, /Pe. are shown on Fig. 6 as the dotted 
and dashed (—— - ——-- ——- - ) lines, with the left hand ends of the 
lines indicating the minimum values of $M _ for which these values hold. 
It will be noted that the computer solution values for all three slender- 
ness ratios are below these limiting values, but appear to approach the 


limits asymptotically as S increases. Once again this is due to the 


20 





bending which has been allowed in the computer solution, but which is 
not considered in the limit analysis. 

The excellent agreement shown between the computer and theoretical 
solutions for the three cases discussed above has been further enhanced 
by three additional checks which have been made on the adequacy of the 
model and computer program to give reliable results. 

In Appendix I, the computer solution for a constant load is com- 
pared with the theoretical solution, given by the secant formula, for the 
same load. The computer solution gives deflections which oscillate about 
some average deflection curve, due to the fact that no damping is in- 
cluded in the model. The "static load" deflections for the computer 
solution have been calculated as the average of the maximum and minimum 
values. It will be seen that the deflection curves compare quite favor- 
ably, even though the computer solution is made with only ten length in- 
crements in the model, instead of the twenty which are used for all failure 
predictions, Also, the dynamic nature of the model response tends to give 
larger deflection values than those predicted by the secant formula. 

Appendix II contains a graphical comparison of the theoretical and 
computer solution for the travel of an elastic strain wave down the 
cclumn. It will be noted that the agreement is very good, in spite of 
the fact that the model used in the computer solution was eccentrically 
loaded and was allowed to bend. The amount of lateral deflection is 
extremely small, however, and does not affect the validity of this check. 

A continuous check, comparing the work done on the column by the 
force pulse with the total energy stored in column as strain (potential) 
energy and kinetic energy, is made during all computer solutions for 


failure values of P, /Pe . These agree within an average value of less 
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than two per cent. The maximum discrepancy found was 6.62 per cent. 

From the data given in Table 1 and the curves of Fig. 6, it appears 
that there are several regions in the failure curves which may be dis- 
tinguished from each other on the basis of the type of failure; that is, 
whether the failure is primarily due to bending strain or to axial strain. 
The type of failure may best be judged from a consideration of the data 
giving centroidal axis strain, €<« , as a fraction of the failure strain, 
0.01. In the case of all three slenderness ratios, beginning with the 
minimum A value, there is an initial region of bending failure, follow- 
ed by a transition region leading, in each case, to the maximum P, /Pe 
values and a region of axial strain failures. This is followed by another 
transition region which leads to the final region of axial strain fail- 
ures, 

Fig. 7 illustrates the lateral deflection change which occurs 
during one of the transitions from bending to axial strain failure. In 
this figure, the lateral deflections at the time of failure have been 
plotted for @ = 2 (€, = 2% of 0.01) and 8B = 7 (€, = 48% of 0.01) for 
S = 150. It will be noted that the maximum lateral deflection for B = 2 
is more than 14 times as large as the maximum for 8 = 7. 

From the data compiled in Appendix V, it is possible to evaluate the 
approximations and assumptions made in developing the equations for the 
system. If this is done, the following is found; 

(a) the maximum error in sin®@; = © is 1.61 per cent 


z 
(b) the maximum error in cos9@j=1- ©; /2 is less than one 
per cent 


(c) the maximum error in l - Add = 1 is less than one per cent 


(d) the maximum error in Av; /)= ©; is 1.56 per cent 
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slenderness ratio 150, for 6 = 2 (P,/P.= 6.60) and 
B= 7 (P,/P= 14.00). 
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(e) 


(£) 


(g) 


i 
9; /2 has a maximum value of 0.05, which is small com- 
pared to 1.00 


the maximum value of |Au;/Z| is 0.05, which is small in 
comparison with 1.00 


the product (e%/2 ) : (Aui/9) is less than 3% of the sum 


(9/2 ) os (Aur/p) 


From this is appears that the assumptions and approximations which have 


been made in deriving the equations are reasonable, and do not produce 


gross errors in the results, 


— OE eagle Age, SE gma 


5, Conclusions 


On the basis of the results discussed in the previous section, the 


following conclusions may be drawn: 


(a) The lumped parameter model in combination with a high speed 


(b) 


(c) 


digital computer provides a powerful tool for the study of 
the dynamic behavior of columns. 

Columns will support loads much greater than the Euler 
load, without failing, when subjected to rapid dynamic 
loading. The maximum load which a column will support is 
dependent upon the type of loading and the slenderness ratio, 
as well as the yield strength of the column material. The 
ability of a column to support large dynamic loads without 
failing is due to the delay in the development of lateral 
deflections caused by lateral inertia. 

Axial inertia effects become significant with this type of 
loading during the first transition from bending failure to 
axial strain failure, and remain important as the rapidity 


of loading is increased, 
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(d) The mode of column failure - that is, whether the 
failure is caused primarily by strain due to bending 
or to axial compression - varies in a distinct manner 


as the rapidity of loading is increased. 
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APPENDIX I 


Comparison of theoretical and computer solution deflections for 
static load, S = 100; e/r = 0.05; P, /P, = 0.85 
The theoretical deflection curve for an eccentrically loaded column 


may be derived from the secant formula as follows: 





\ 
Umax . € Bu \2 
—* EL sec (BS -1 
where a is the maximum deflection from the original position. Now, 
P, = 0.85P. = 0.85 Tt EI/L* , so that 
VV 
Vmax . ©€& (ee 
pax = £ [sec ( 9 ‘ 


However, e/L = e/(rS), and e/r = 0.05, S = 100. 





Therefore 
Umax — 0-05 | . | 
cL 105 7.2055 


In order to get vo /L on a percentage basis, multiply both sides by 100, 
giving 

mtnecht 2) = 0.05 [ 7.2055] = 0,360275 
The remainder of the theoretical deflection curve may now be calculated 
from 

v/L(%) = 100(v’_ | /L) sin( T% x’ /L’) - 100 (e/L) 


which is derived from a consideration of the sketch below. 
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The theoretical deflections computed from this equation and the 
computer solution deflections are tabulated in the accompanying table 
(Table 2), and are compared graphically in Fig. 8. The computer solu- 
tion does not yield a constant set of deflections, but rather, gives 
values which oscillate about some average deflection curve. This is 
due to the fact that no damping is included in the mathematical model. 
The static load deflections have been computed as the average of the 


maximum and minimum values given by the computer solution. 
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Table 2 
Static Load Comparison 


S = 100; e/r = 0.05; P, /P, = 0.85 








a 
0.85 
Owes. 
SS 
ie 
0 
1 Za 3 oa 5 6 
Time (T) 
Computer Solution Force-Time Historv 
Per Cent Theoretical Maximum Minimum Average 
Length Lateral Defl. Computed Computed Computed 
(x/L x 100) eri) Lateral Defl. |Lateral Defl. | Lateral Defl. 


(% L) (% L) (% L) 


) C 


10 Oe 12 0.209 Goo 1 On 20 
20 o.213 0.396 O57 0.226 


30 0.292 0. 540 0.076 0. 308 
40 0. 342 0.635 0.086 0.361 
50 0. 360 0.670 0.092 0.381 


60 0.635 0.086 0. 361 

70 0.076 0. 308 

80 pa2ig 0.057 0.226 

90 0.112 0.031 0.120 
0 


100 | 0 
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APPENDIX II 
Comparison of theoretical and computer sclutions for an elastic 
strain wave in the column. S = 50, 4 = 4, A, = 10, T = 0.0392699. 
The theoretical solution for the strain at any point in the column 
may be derived as follows: 
At a time T, the force at the loaded end of the column is 
A. = S = A, sin (1tgT) 


and, neglecting the effects of lateral response motion, the force at some 


distance x from the loaded end is 


A: = = = A, sin [8 (T - oa) 


where x/s is the (real) time required for an elastic wave to travel a 


distance x. 


Now, since F = EA/S’ , 


P. = (=) Ag Sin [te (t- = ) | 


{ 


However, € = P/EA, so that 
As c: i 
Ee; = ; sin |r = 3 
L Sg B ( 7 
Since @, = 2LS/tt s, this may be written as 
= ome [ | 
— Sie ne T= 
=> = A(t aS) 
which, upon substitution of A = 10, S = 50; B = 4, and T = 0.0392699, be- 











comes 
€, = 4x10> sin | 416 (0.0342694 - a 
The theoretical strain values have been computed from this equa- 
tion and are compared graphically with the values from the computer solu- 
tion in Figs 9. 


It will be recognized in the last equation for €i that x must reach 


a value of 1.25L before €:; = 0. This means that the pulse has reached the 
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pinned end of the column and has been reflected back a distance of 0.25L. 
Thus, the values of €; computed for 1 <= x/L —1.25 must be added to the 
values calculated for 1=>x/L = 0.75 to give the resultant strain in the 


region 0.75 <= x/L —1. 
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APPENDIX III 


Computer program notation 


Computer 
Program 
Term 
ALPHA (TI) 
ANGL (I) 
BENDM (TI) 
BETA (I) 
DDX (I) 
DDY (TI) 
ECC 
EXTSTR (I) 
FORCE 
FAC] 

FACZ 


FAC 3 


Corresponding 
symbol in the 


basic eqs. 

Ay 

oy 

Yn: 

Bi 
d*(Yi/y) aT? 
d*(¥/0)/dT* 
I 


e 


rp 
Ey (max) 
A, 


x 


eS 
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Computer Corresponding 
program symbol in the 
term basic eqs. 
INTV n 

J n+ 1 

RHO p 

SLRT S 
STRAIN (I) Cnr 

TIME T 
TIMEINT AT 
xacch (1) (4¢Yp) fat? 
XDEFL (I) (Ui/Q) 
XVEL (1) d (4i/9) [at 
yacch (1) = & (“*/2) Jat? 
YDEFL (I) (vi /2) 
YVEL (I) A(Y /g) aT 





PROGRAM COL/ 


PROGRAM TO OETERMINF COLUMN -BEHAVIOR WHEN HALF SINE PULSE LOAD ILS 
APPLIED WITH SMALL ECCENTRICITY 10 INITIALLY STRAIGHT COLUMN 


100 
110 


190 


ODIMENSION ALPHA(60),BFT4 (60), XACCL(60),YACCL(60), XVEL(60), 
TYVEL(60)_9 XDEFL( 60), YDEFL(60),BENDM(60) ,ANGL(60)5STRAIN( 60), 
2 ull), P(l)s EXISTRE4Q) 
O COMMON INTV,Je TIMEINFE,XACCL XVEL + XOEFL,VACCL,YVEL,YDEFL, 

| AINTV,RHO,ANGL,BENDM,ALPHA,BETA,TIME,ECC,FAC2,FORCE,C1,02,C5;, 
2 WORK | 
INTV = 20 

READ 100, SLRT-, FACI, FAC2, FACS, FORCE, ECC 

READ 110, A 

FORMAT (569.2, §20.6) 

FORMAT (€9.2) 

J = INTV + 1 

AAINTV = INTV 

RHO = SLRT/AINTYV 


DO 200 I = 1,J 
ANGL(I) = O. 
XVEL(I) = O. 


YVEL(I) = QO. 
XDEFL(I) =0.- 
YDEFL(I) = Q. 


200° BENOM(]), = 0. 
- FORCE = FORCE - A * 9.87 
WORK = 0. 
TIME = 0. 


210 
220 


2300ALPHA(T) 


ALPHA( 1) = QO. 
FAC} = 0.40 * FaAC2 
TIMEINT = (1.570796) /S((AINTV8#2) #RHOFFAC I) 


Cl = (0.405285) #(AINTV##2) 
C2 = 2. * Cl 
C3 = TIMEINT/2-. 


DO 220 I = 2,INTV 

BENDM( I) = (AINSI V##2) #(ANGLO(ITY -ANGLOII+1)) 

DO 230 I 29J 

~(AINTV##2)#(RHO##2)*%(XDOEFL( IT) ~XDEFL(I-1) 
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235 


240. 


250 
255 


256 


1 +(ANGLOI)##2)/2.) 


BENOM(J) = ALPHA(J) #*°ECC 

Dop2somt = 2,0 

BETA(I) = ANGL(1L) © ALPHACT) - (BENDM(I) -BENDM(I-1)) 
BETA( 1) = BETA(2) 

ALPHA(J+1) = ALPHACd) 

BETA(J+1) = BETA(J) 


DO 240 I = 2,J 
XACCL(I) = CIl*#(ALPHA( I) - ALPHA(I+4+1)) 
YACCL(I) = Cl * (BETACI) ~-BETACI+1)) 
XACCL( 1) = C2 # (ALPHAC1) -ALPHA(2)) 
YACCL(1) = O. 

CALL SIGMA 

DO 255 I = 23J 

STRAIN(I) = -CXDEFL(I) -~XDEFLCI-1V) +(ANGL(T) #*#2)/2.) 
STRAIN( 1) = OQ. 

DO 256 I = 2,INTV 

BENOM(IT) = CAINTV#¥#2) *#(ANGL(ITI-ANGL(I+1)) 


DO 260 I = 2,J 


Z2600EXTSTRII)= ABSF ( STRAINGI)) + € Of fS/(RHO#AINTV #42) ) 


900 


950 


905 
910 


915 


* ABSF (BENDM(1) + BENDM(I-1)) 
EXTSTR( 1) = O.- 
I= 1 
I=I+ | 
IF ( I - J ) 905,905,950 
IF( TIME -— I1-/FAC3 ) 210,710,710 
IF ¢€ EXTSTRCI) — 0-01 ) 900,910,910 
SUMKEN = 0. 


VEN = 0. 

_BENEN = ©. 
COMEN = 0. 
TOTEN = OG. 


DO 915 1 = 2yINTV 

VEN=VEN#( 1223370 1/AINTV# #2) ® (XVEL(T) ##2+YVEL(T) #82) 
BENEN=(0.50/AINTV##2)*(BENDM(1)##2) + BENEN 

SUMKEN = VEN + (0.616850/AINTV##2) #(XVEL( 1) *#2) 

DO 916 I = 2eJ 
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916 COMEN=COMEN+ (0.50) #(AINTV##2)#(RHO##2) #(STRAIN( 1) a2) 
TOTEN = SUMKEN + BENEN + COMEN 
PRINT 300 
Z000FORMAT (6H) SLRT(X,5HFORCE?X, 3HECCTX gp 4HEAC26X%y (HTIMEINT3Xy 
IUHINTV//) , 
PRINT 400 (SLRT, FORGE, EGG, FAG2, TIMEENTs ENTV) 
400 FORMAT ( IPE9.2, IP4UE1I.2, IT////) 
PRINT 500 
500 FORMAT (20H TIME//) 
PRINT 510, (TIME) 
510 FORMAT (1PE20.7, 1H //) 
PRINT 550 — 
550 FORMAT (5H ISX, 5HXDEFL 8X ,SHYDEFLTX,6HEXTSTR//) 
PRINT 600, { I» XDEFL(I), YOEFL(I), EXTSTR(I), I = yd) 
600 FORMAT (15, 1P3613.2) 
PRINT 601 
601 FORMAT ( HO ) 
PRINT 605 
605 FORMAT (20H TOTENI6X,4HWORK//) 
PRINT 606 ( TOTEN, WORK ) 
606 FORMAT ( 1P2E20.7 ) 
PRINT 800 
800 FORMAT (1H1)- 
GO TO 190 
710 STOP 
END 
SUBROUTINE SIGMA 
ODIMENSION ALPHA(60),RETA(60), XACCL(60),YACCL(60), XVEL(60), 
1YVEL (60), XDEFL( 60), YDEFL(60),BENOM( 60) ,ANGL(60),STRAIN( 60), 
2U(1),P4 1), ODX(60), DDY(60), XSPD{(60), YSPD(60), X(60), Y(60) 
O COMMON INTV,J,TIMEINT,XACCL,XVEL,XDEFL», YVACCL,YVEL,YDEFL, 
1 AINTV,RHO,ANGL »yBENDM,AL PHA, BETA, TIME, ECC, FAC2,FORCE,CI,C2,C3; 
2 WORK 
UCT) = XDEFL(1) 
P( 1) = ALPHA(1) 
00 10 I = IyJ 
XSPD(T) XVEL{I) +# TIMEINT * XACCL(I) 


tl 
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20 


25 


26 


27 
28 


30 


40 


~~ 


50 


55 


YSPO(I) = YVELCI) + TEMEINT * YACCL(I) 


X(I) = XOEFL(T) #C3 ® (XVELCI) + XSPDCI)) 
YCI) = YDEFL(I) + C3 * (YVEL(I) + YSPD(I)) 

DO 20 I = 2sJ 

ANGL(T) = YCI) - YCI - 1) 


TIME = TIME + TIMEINT 
IF (TIME #® FAC2 <1.) 253 263,276 


ALPHA( 1) = FORCE * SINF(FAC2 * TIME * 3.141593) 

GO TO 27 

ALPHA( 1) = O. 

DO 28 I = 2,J 

ALPHA(I) = -(AINTV##2)*(RHO##2)#(X(1L)-XCI-1T)+0ANGL( I) ##2)/2.) 


ALPHA(J+1) = ALPHA(J) 

DO 30 I = 2,1INTV 

BENDM(I) = CAINTV##2) 4(ANGL(OI) -ANGLOUI+1)) 

BENDM( 1) = ALPHA( 1) * ECC 

BENDM( J) ALPHA(J) # ECC . 

DO 4O I = 2,J 

BETACI) = (ANGLO I)*ALPHA(I)) - (BENDM(IT; - BENDM(I-1)) 
BETAC]) = BETA(2) 

BETA(J+1) = BETA(J) 

ODX( 1) 
ODY( 1) = O.- 
DO 50 I= 2.J 


C2#(ALPHA( 1) - ALPHA(2)) 


DDX(1) = Cl®(ALPHA(T) -ALPHA(T+1)) 
ODY(I) = Cle(BETA(I) -BETACI+1)) 

pO 55 I= Ip 

XSPD(I1) = XVEL(I) + C3*(XACCL(I) + DDX(I)) 
YSPD(I) = YVEL(L) + C3*(YACCLII) + DDY(1)) 
XDEFL(I) = XDEFL(I) + C3e(XVEL(I) + XSPDC1)) 


YDOEFL(I) = YUDEFL(I) + CS3#(YVEL(I) +YSPD(T)) 


XVEL(I) = XSPD(1) 
YVEL(I) = YSPD(1) 
mpor60) 1 = 2,J 
ANGL(I) = YOEFL(I) - YDEFL(I-1) 


60 


WORK = WORK + ( XDEFL(1) -UCT))*#( PC 1) + ALPHA(1))/2. 
RETURN 
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APPENDIX IV 


Failure Data 


et ee 


Tabulated in the following pages are the failure data for the 
three slenderness ratios for all values of @ which were investigated, 
These data are the x and y direction deflections and the extreme fiber 
strains at the time of failure, for the minimum failure values of P, /P.. 
The strain values given are the strains half-way between the point for 
which they are listed and the previous point. Thus, a strain listed 
opposite 50 per cent of the length is actually the strain at 47.5 per 


cent of the length. 
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Failure Deflection and Strain Data for Minimum Failure 


Value jof P, /P, 
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Bare Deflection and Strain Data for Minimwn Fe-ilure 
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Failure Deflection and Strain Data for Minimum Failure 


Value joe P,/Pe 
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Failure Deflection and Strain Data for Minimum Failure 
Value of P, oe 
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Failure Defleotion and Strain Data for Minimum Failure 


Value of P, /P. 
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Failure Deflection and Strain Data for Minimum eens 


Value of P, /P, 
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Failure Deflection and Strain Data for Minimum Failure 


Value of P, /P. 
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Failure Deflection and Strain Data for Minimum Failure 


Value Oe Pi ve 
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Failure Deflection and Strain Data for Minimun Failure 
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Failure Deflection and Strain Deta for Minimum Fa:ilure 


Value jof P, /P- 
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Failure Deflection and Strain Data for Minimum Failure 


Value of P, /P, 
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Failure Deflection and Strain Date for Minimum Bajilure 
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Failure Deflection and Strain Data for Minimum Failure 


Value of P,/P. 
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Tailure Deflection and Strain Data for Minimum Failure 


Value of P, /P, 
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Failure Deflection and Strain Data for Minimum Failure 
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Failure Deflection and Strain Deta for Minimum Failure 
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Failure Deflection and Strain Data for Minimum Failure 
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Failure Deflection and Strain Deta for Minimum Failure 


Value of P,/P, 


5 = /00 1? = (4.478 
AT = 9 82 vjo77 
n =Zo 


Lateral Eoctreme 
deflection Yiber 
(vi /0) .  gtrain 





a ee ee 
Be 93 

8. 83 

8.67 

8 hb 

8.19 
7.87 
7.50 
7.09 
irr: 

be / 
wl 
P. wo, OT Ee 
$60 1) 
_whJoG WIM) 
3.36 
80] le CT 
| 


ol, OF 
Te: ee ee Ne De ey ee. TT ET Rae hl weIay Qs 






7 /e +5 X /o7* 


90 _h $0 r 
95 Oe FIR ft -&. kf Goo 


- 100 O O "3G 4 


59 


- 
PS. 





: SS 


[45 | 6 es 


Failure Deflection and Strain Data for Minimum Failure 
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Failure Deflection and Strain Data for Minimum Failure 
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Failure Deflection and Strain Data for Minimwn Failure 
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Failure Deflection and Strain Data for Minimum Failure 
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Failure Deflection and Strain Data for Minimum Failure 
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Failure Deflection and Strain Data for Minimum Failure 
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Failure Deflection and Strain Data for Minimun Failure 
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Failure Deflection and Strain Data for Minimum Failure 


Value of P, /P. 
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and Strain Date for Minimum Failure 
Value of P, /P, 
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Failure Deflection and Strain Data for Minimum Failure 


Value of P, /P. 
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Failure Deflection and Strain Date for Minimum Failure 
Value of P, /Pe 
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_ Failure Deflection and Strain Data for Minimum Failure 







Value of P, /P. 
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Failure Deflection and Strain Data for Minimum Failure 


Value of P,/P. 
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Failure Deflection and Strain Data for Minimum Failure 


Value of P, [Pe 
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Failure Deflection and Strain Data for Minimum Failure 


Walue Sof 2 7 Pe 
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Failure Deflection and Strain Data for Minimum Failure 


Value ‘of P, /P. 
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